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Quantum information processing has greatly increased interest in the phenomenon of 
environmentally-induced decoherence. The spin boson model is widely used to study the interaction 
between a spin-modelling a quantum particle moving in a double well potential-and its environment- 
modelled by a heat bath of harmonic oscillators. This paper extends a previous analysis of the static 
spin boson study to the driven spin boson case, with the derivation of an exact integro-differential 
equation for the time evolution of the propagator of the reduced spin density matrix. This is the first 
main result. By specializing to weak damping we then obtain the next result, a set of Bloch-Redheld 
equations for the equilibrium fixed spin initial condition. Finally we show that these equations can 
be used to solve the classic dissipative Landau-Zener problem and illustrate these solutions for the 
weak damping case. The effect of dissipation is seen to be minimised as the speed of passage is 
increased, implying that qubits need to be switched as fast as possible. 



-a 
c 

o 
o 



> 
^1- 



o 

-I— > 



I 

C 

O 

o 



X 



I. INTRODUCTION 

A new era of quantum information processing experi- 
ments is beginning (e.g. 0,0)01) with the consequence 
that battling decoherence (the destruction by environ- 
mental degrees of freedom of the phase coherence be- 
tween superposed quantum states) has become a crucial 
task for designing the fault-tolerant quantum computer 
The engineering significance now attaching to de- 
coherence gives new urgency to an already fundamental 
theoretical question (e.g. S 0) : wna t is the effect 
of an environmental heat bath on a quantum two-level 
system (qubit)? 

As well as its importance for the experimental realiza- 
tion of a qubit [aSEy]' the current interest in this prob- 
lem is due to at least two other factors. One is its gen- 
eral importance as a testbed for com pet ing conceptions of 
dissipative quantum mechanics [ll], [l^] • A second is the 
bridge it offers between quantum and classical stochastic 
systems, for example as a quantum rather than classical 
two level stochastic resonance system (e.g. 0,^3)- 

A particularly well-studied framework for the problem 
is the spin-boson Hamiltonian (e.g. 0,01) i n which the 
two level system is modelled by a spin, the environmental 
heat bath by quantum harmonic oscillators, and the spin 
is coupled to each bath oscillator independently. This 
Hamiltonian is easily adaptable to the driven spin-boson 
case where the two level system is subjected to an ex- 
ternal force. In the strong coupling and high tempera- 
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ture bath limits much progress has been made, especially 
for the case of harmonic driving which defines the quan- 
tum stochastic resonance problem 0, 0] . Hartmann et 
al. 0| have, however, recently emphasized that weak 
coupling and low temperatures remain much less well ex- 
plored, and that although formal solutions for the result- 
ing spin dynamics do now exist, they still in practise have 
to be used perturbatively. 

Hartmann et al. |14j sought to make progress by point- 
ing out similarities of the driven spin boson problem to 
counterparts in quantum optics and solid state physics 
where Bloch-Redfield equations are used. They exhibited 
a set of Bloch-Redfield equations (their Eqn. 4) obtained 
by projection operator methods. 

In this Paper we show that an earlier path integral- 
based derivation of a similar set of equations (la ] for 
the time independent spin-boson Hamiltonian can be ex- 
tended to the driven case (see also |l6|). This allows us 
to derive (section [HJ an exact integro-differential evolu- 
tion equation for the propagator of the reduced density 
matrix, and is the first of of three principal results of the 
Paper. We then (section IIII(I make the weak coupling 
assumption to give an independent derivation of a set of 
Bloch-Redfield equations for the driven spin boson case; 
the second main result of the Paper. We note in pass- 
ing that this refutes the assertion by Zhang that the 
methodology used in |l5j was unsuitable for time depen- 
dent problems. 

We choose an initial condition in which the environ- 
ment equilibrates about a fixed spin, which is subse- 
quently released at time to- We find, in consequence, a 
difference from the Bloch-Redfield equations of 0], but 
only in the terms arising from dissipation. This difference 
disappears if we follow previous authors in assuming spin 
and bath to be completely uncoupled at to- In either case 
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our terms due to fluctuation are identical with those seen 
by 0. 

The main application we study for our Bloch-Redfield 
equations is the dissipative version of the classic Landau- 
Zener-Stuckelberg fLZS) nonadiabatic level crossing 
problem (e.g. |a 111 Ha l20l| ). This describes the situ- 
ation where the energy levels of a quantum two-level sys- 
tem are brought together over time by an external force, 
enabling transitions to take place. Examples of this ubiq- 
uitous problem are found in atomic collisions, chemical 
reaction kinetics, biophysics, solar neutrino oscillations 
and NMR (see Ao and Rammer for references and 
additional examples). In the present case we will take 
the LZS Hamiltonian with a fixed tunnelling term and 
linearly time- varying bias, and study the effects of an en- 
vironment consisting of harmonic oscillators. The initial 
study of of the dissipative LZS problem by Ao and Ram- 
mer |18| has since been supplemented by other analytic 
and computational methods (e.g. |a|20|,|2l|]). In one case 
|2l| these benefit from the recent feasibility of computa- 
tionally following ss 120000 basis states; but they have 
all given results for the influence of the bath, for inter- 
mediate speeds of passage, which disagree with those of 
[l8| . Here HHIf) we specialize the Bloch-Redfield equa- 
tions we have derived to the LZS problem in the ohmic 
weak damping case, which greatly simplifies the dissipa- 
tion term (as also found by These equations for the 
dissipative LZS problem may aid study of topical issues 
such as the merits of different models of the noise power 
spectral density J{uo) of the heat bath || or of alterna- 
tive functional forms for the time dependent bias in the 
LZS problem (c.f. Chapter 6 of 

By scaling the Bloch-Redfield equations to the char- 
acteristic tunnelling time, we show that dissipation is 
significant even at zero temperature, for finite speeds of 
passage, agreeing with newer authors, but that its im- 
pact is inversely proportional to the passage speed. This 
implies the intuitively reasonable conclusion that qubits 
need to be switched as fast as possible to minimise dis- 
sipation effects. We illustrate our findings numerically, 
in the zero temperature case where the fluctuating force 
terms greatly simplify in addition to the more general 
ohmic simplification of the dissipation term. Damping is 
clearly seen, reminiscent of that seen for the high tem- 
perature case in 0] 



II. EXACT INTEGRO-DIFFERENTIAL 
EVOLUTION EQUATION FOR SPIN 
PROPAGATOR OF DRIVEN SPIN-BOSON 
HAMILTONIAN WITH EQUILIBRIUM FIXED 
SPIN INITIAL CONDITION 

First, in subsection III Al we describe the driven spin- 
boson system, definin g th e Hamiltonian and our nota- 
tion (based on [El ITEl l26|L In subsection III Bl we write 
down the reduced spin density matrix as a path inte- 
gral, to which (subsection III Cfl we apply the equilibrium 



fixed spin initial condition. We can then (subsection lll Dl 
and lll Kfl integrate out the bath variables to find the effect 
of the bath on the spin by means of the Feyman- Vernon 
influence functional for this initial condition. This en- 
ables us to generalize the results of Waxman |l5j to the 
driven spin boson case, arriving at our first main result 
( 1471) , a general time evolution integro-differential equa- 
tion for the propagator K of the reduced density matrix 

In this section we make only two main assumptions 
in addition to those of the spin-boson model itself: i) 
the widely-employed (e.g. [l8|]) choice of a factorizing 
initial density matrix and ii) use of an "equilibrium fixed 
spin" initial condition where the spin is assumed to be 
fixed "up" and in equilibrium with the bath at to- To aid 
subsequent comparison with Hartmann et al. |14| we also 
give the equivalent influence function for the case when 
spin and bath are uncoupled at to. 

A. The Driven Spin Boson problem 

The driven spin-boson Hamiltonian H is: 

H = H S + H I + H B (1) 

where 

Hs = -—o x + -^a z (2) 

where o~i with i = x,y, z are Pauli spin operators; —A and 
e(t) are angular frequencies corresponding to off-diagonal 
(tunnelling) and on-diagonal (bias) spin matrix elements 
respectively; and the heat bath is represented by a set of 
harmonic oscillators of mass m a , angular frequency u> a , 
momentum p a and position coordinate x a . The oscilla- 
tors are coupled independently to the spin co-ordinate 
with strength measured by the set {c Q } while qo mea- 
sures the distance between the left and right potential 
wells. 

In contrast to most previous work (e.g. |15^ we will 
use an initial condition where at times — oo to to the spin 
has been held fixed and the heat bath has equilibrated 
around it; a choice which requires explanation. For rea- 
sons of mathematical simplicity we want a factorizing 
initial full density matrix. Given this, we have selected 
what we believe to be the boundary condition that avoids 
unnatural transients associated with the time evolution 
of the complete system after to . Any other choice for the 
bath's initial density matrix would be expected to lead 
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to transients because at time to both the bath and the 
spin would have dynamics driven by the interaction be- 
tween them. Such an initial condition can be viewed as 
the result of having e(t) very large and negative for times 
t < tQ. 

We thus take 



H(t = to) = 



E 



2m a 



m a u)l 



qo c 



(5) 

which incorporates the initial condition for the combined 
interaction and bath Hamiltonian at time to- This results 
from (JIJ with the eigenvalue of a z set at +1. The extra 
term relative to namely 



go c 



2vn a Lol 



cancels with its equivalent in the partition function (the 
denominator of l|15(l . 



B. The density matrix as a path integral 

We follow the method of 15] but the driven spin- 
boson Hamiltonian (JTJ replaces the static spin-boson 
Hamiltonian given as his equation (1). We start by ex- 
pressing the density matrix of the combined system as 
a double path integral. The density operator p obeys 
ihdp/dt — [H,p], so its time evolution is given by 
p(t) = U(t,to)p(to)U~ 1 (t,to) with U a unitary time evo- 
lution operator. 

Hence with \xa > an eigenstate of the oscillators' co- 
ordinate operators x and the spin operator a z we have 



< x 1 <j 1 \p(t)\x 2 (J2 > = p{xia 1 ,x 2 <T2;t) 

= EE j dXs j dXA < Xl<Tl I U ( t ' t o) \ X 3 a 3 > 



<J3 (7 4 

r-l 



x < x 3 a 3 \p(t ) \x 4 a 4 >< x 4 a 4 \U (t,to)\x 2 a- 2 > (6) 



where we have inserted two complete sets of states and 
boldface type distinguishes a state of all N oscillators 
(x = {x a }). We now assume that the initial density ma- 
trix factors into an oscillator dependent part and a spin 
dependent part. We are only interested in the behaviour 
of the spin so work with the reduced density matrix, ob- 
tained by integrating over the bath variables. We then 
define the reduced density matrix p by tracing over the 
bath variables, and normalizing to the free oscillator par- 
tition function. We have 



p(aia 2 t) = ^ ^ K{aia 2 t\(T^(JAto)p{(T Z (j A to) 



(7) 



with 



S[<j,x] = S s [(t] +Si[a,x] +S B [x] (10) 
where Ss,Si,Sb are the actions corresponding to the 
spin, interaction and bath Hamiltonian operators Hs, Hi 
and Hb respectively. The notation 



Xit 



d[x] 



(11) 



Our initial problem is to determine the propagator K: ,i u ,i , • • , _ , , . , 

1 * f & means the sum over all paths beginning at 033 at time to 

/[■ r and ending at X\ at time t, and the spin path integral is 

dx' / dx 3 / dx4(aj / <7i|f7(t,*o)|sc3<7'3) explained more fully below. 

x p(x3X4to)(x4<j4\U~ 1 (t,to)\x' (T2) (8) Combining the above with a similar expression for the 

backward propagator and with the initial bath density 
the effective time evolution operator for the reduced den- matrix we obtain 
sity matrix. 

We now write the matrix element of the time evolu- 
tion operator U(t, to) (the forward propagator) as a path 
integral 



(xiai\U(t,t )\x 3 a 3 ) = / d[x] / d[a] exp I -S[tr, x] 

JX 3 to Ja 3 t V" 

(9) 



/■oil /"«2l 

K{aia 2 t\az<T4to) = d[a] d[v}A n [a]A*Ju] F[a, v] 

(12) 



where 
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F[a, v] = 



/p pX' t pX't 

dx 3 / dx± / d[y] / d[x] 
J <J X 4 £q J X3 t f j 

xe HSi[^]-S 1 [^y]+Sn[x]~S B [y]) 

I 



(13) 



is the influence functional containing all effects of the description of the oscillator bath at time to- We assume 
bath on the spin system. thermal equilibrium so 



d[a] exp — S[cr 
n 



dt r , 



*2 



with the sum[[,]] running over all numbers of flips (points 
at which <r goes from —1 to +1 or vice versa) consistent 
with the initial conditions, defines both the measure d[a] 
and the amplitude A n [a] associated with a path [a] con- 
taining n spin flips. 



C. Application of equilibrium fixed spin initial 
condition 



p(x 3 x 4 t ) =< x 3 \ 



(15) 



in which Hq denotes H(t — to). The partition function 
is given by Z Q = Tre-P Ho and the inverse temperature 

by/3: 



1 

k B T- 



We now calculate the influence functional and so have 
to specify p(x 3 X4 : to), which corresponds physically to the 



We find that the path integral part of the influence 
functional (|13f) factors into a product of terms like 



J 



d[y) / d[x]<x 3 \ — |x 4 > exp~(S a i[a,x] - S a j[v,y] + S aB [x] - S aB [ 

Jxztr, z n 



X4,t 



(16) 



r 



one for each [[]] oscillator, x 3 and X4 are the initial states we first evaluate the forward propagator for the bath vari- 
for the forward and backward paths. Noting that the ables, a standard result [23| : 
(single particle) bath and interaction actions are 



S a i[a, x] +S aB [x] 



1 , ( m a x 2 a m a u 2 a x 2 a c a x a q a 
du 



In 



2 
(17) 



J 



x'i 



X3to 



d[x] exp - (S a i [a,x}+ S aB [x] 



2nih sin u> a (t — to) 
im a uj a 



x exp 



21% sin u) a (t — to) 



-B 



(18) 



with 



-x'f a (t - t ) + X 3 ga(t - t ) 



B = (x\ + (x'f) cosw Q (t - t ) - 22/23 
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o o f~~) I du I dv a (u)<t(v) sin ui a (t — u) sin uj a (v — tcf^u — — Ca< ^° / g?m <r(u) sinu a (i — u). (20) 



f a = / g?w <t(m) sin w Q (w — to), (19) propagator gives 

J to 



Combining this result with the its equivalent backward 



x't r-x't 

d[x] / d[y]exp 

x a t J X4,t 



■ (S a I [cr, x] - S a l [v, y] + S a B [x] - S aB [ 



■ exp 



2ttH sinui a (t — to) \2hsinuj a (t — to) 



-C 



(21) 



where 



C = (x\ — x\) cosuo a (t — to) — 2x'(x 3 — Xi) 

-x'[fa(t - t ) - f v (t - t )} + x 3 g a (t - t ) - Xig v (t - t Q ) 

-44 f du r dv 

miu>a Jt J t 

x [a(u)a(v) — v(u)v(v)\ sin tu a (t — u) sin uj a (v — to) (22) 

I 



and the functions /„ , <?„ are f„ , g a but with v{u) replac- 
ing a(u). All of the above is still just for one oscillator 
hence the use of the a index. We now do the integral over 



x' in equation f!13[). which yields a delta function, so the 
path integral part of the influence functional becomes 



X t rX t • 

dx' I d[x] d[y]exp-(S a i[a, x] - S aI [v,y] + S aB [x] - S aB [y] 

x 3 t Jx 4 t " 



X6 



X 



{2(x 3 - x 4 ) - (/„ - /„)} 



, iX 
exp ( — D 



V 



(23) 



where 



and 



X = 



Ti sin uja{t — to) 



(24) 



D = (xl~ xl) cos u> a (t- to) + x 3 ga-(t- t Q ) - Xig u (t- t ) 



2 2 



du / dv 

to 



x [a(u)a(v) — v(u)u(v)] sin w a (t — u) sinuj a (v — to)- 
I 



(25) 



We note that the above expressions depend on the sign of the coupling constant c a via the j aM and g a v terms, but 
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not the initial value of a (because (|23|l is a product of two 
propagators). The initial value we chose for a, (= +q /2) 
did, however, appear in the thermalized density matrix 
(US)) for the bath. 

We recall from ( I13|) and ( I15fl the full influence func- 
tional including all the oscillators. The path integrals 

factor similarly, so JxJ° d[x] is a product of single par- 
ticle path integrals over all paths that the individual x a 
can follow consistent with the set {x a } being x' at time t, 
and £c 3 at t . If we can "fold" in the product of the traces 
of single-particle propagators given by f2"3"|) and f2"5|) we 
can then multiply by p(x 3 x^to) and do the integrals over 
x% and X4 to obtain the influence functional. 



D. Evaluation of the influence functional in 
equilibrium fixed spin case 

We want to evaluate a non-standard influence func- 
tional, corresponding to the equilibrium fixed spin initial 
condition. We do so by replacing it by the influence func- 
tional for a system with the spin and bath uncoupled at 
to multiplied by a phase due to the fixed spin boundary 
condition. To derive this result we now go to sum and 
difference coordinates z\ — X $+ X4 - and z 2 = x 3 — X4. 

If we define F a [a, v\ , the single oscillator influence func- 
tional, by F[a, v] = Yl a F a [<^,^], then using (jTTj|l . and 
H a o of the form given in (JSJ, we have 



f ( 1 \ z 2 e'f 311 " z 2 

FJpM = J dzidz 2 S[zi- -(f a - f v )\ < zi + —\ — - — \z\-—> 

( im a uj a \ 

x exp — — - -M (26) 

y \2hsmu) a (t-t ) J K ' 

I 

and x [a(u)a(v) — v(u)v(v)] sinw Q (t — u) snxu) a (v — toJ%7) 

M = z\(f a - f v ) cosuj a (t - t ) + zi(g a - g v ) + y(&r + 9v) 

n 2 2 ft r-u where the delta function has been used to replace x\ — x\ 



du dv by Zi(f (7 — fv)- We now consider 

J to J to 



'■a UJ a J t J t 
I 



z 2 1 exp -f3H aQ z 2 Posc 
< zi + —I - \zi - — >= -=-(z\ + z 2 /2 + A, Zl - z 2 /2 + A) (28) 

l ZjO £ ^0 



where If we now write z\ + A = z 1 we get 

A = qpc a = qoc a a(t ) ^ 
2m a u 2 a 2m a uj 1 l 

the amount by which the coordinate x a has to be dis- F ^ ^ = FaA=Q ^ „] exp (p) ( 31 ) 

placed to go from the free single oscillator density matrix 
p osc to the fixed spin case. In fact p osc is just exp (—(3H a ) 
with: 

v 2 1 

H « = 7^- + T^^lxl. (30) with 



p = 2nZ7-to) { - A) [(/ * - cos ^ (< - <o) + ^ ~ ^ ' (32) 



and F a A=o the conventional influence function (|37|l . We also have that 

(fa - fv) COSUJ a (t - t ) + (g a - Qv) 
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c a qo 



du[a(u) — smu a {t — to) cosw Q (u — to) 



so the full influence functional is 



and 



F[<J,v] = f 7J F aA = ) exp ° ° A I du[a(u) - cos iu a (u - t ) 



2h 



i<lo 



= Fa=o exp — a (to) j du[a(u) - v(u)]2Qi(u - to) 



(33) 



r 



where 



Fa=o — 77 F a A=o, 



is the bath spectral function describing the oscillator heat 
bath. Equation l|33|l should be compared with the ex- 

(34) 

pression arrived at by Leggett et a/|5J] as their equations 
(B.1.9) and (B.1.10). All we need now is the standard 
result 



l( u ) = ^ o C " a COSUJ a U 
^ 2m a u>^ 



dcu J(u>) 



■ COSLUU. 



IT LU 



(35) 



J(uu) 7T^ C 2 a 



2 ^ m a u>%. 



6(u> — U> a ) 



(36) 



Fa=o[o", = exp-— j du j dv[a(u) — u(u)][a(v) +v(v)]Q x (u— v) 



e 



x exp — 57 du dv[a(u) — v(u)][a(v) — v(v)]Q2(u — v) (37) 

I 

where t = f3h. We also have the conventionally defined 

■ ■ ' Lu olll Lull I O J I 



correlation function Q2 for the fluctuating force (see also J Q it lu 

Appendix[Bj and the retarded resistance function Q x (see 
0) given by 

Qi(u) = / — J(w)coth — coswm, (38) 
Jo t 2 

and 

7 respectively. The propagator for the reduced spin density 

Qi(u) = —Qi(u) matrix in ( 0) becomes 

du 



pcrit rait 

K(a 1 a 2 t\a^a A to) = J d[a] / d[v]A n [a]A* m [v\F[<x,v] 

C<Tit poit 

d[a] / d[v]A n [a]A* m [v]F A=0 [<r, v] 
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x exp^cr(t ; 



duQi(u — to)[a(u) — v{u)]. 



(40) 



to 



We see that the equilibrium fixed spin initial condition by the spin boundary condition. From (|37fl and i' ; >.'>ll we 
has altered only the imaginary (dissipative) first factor in have: 
equation 7 of |l5j , Waxman's expression for the influence 
functional. The second (real) factor due to fluctuations 
is unchanged. 



E. Simplification of the fixed spin influence 
functional 

We can simplify the first part of the influence func- 
tional considerably by using the extra factor introduced 



F[a,v] = exp{-a(t )(y) 2 y t du[cr(u) - i/(«)]2Q x (« - 1 )} 
iq 2 f* f u 

xexp{--^- du dv[a(u) - v{u)][a(v) + v(v)]Q 1 (u - v)} 
q 2 [* f u 

x exp 2. du dv\a(u) — v(u)\ \a(v) — v(v)]Q2(u — v). 

Ah Jt J t 



(41) 



r 



We note that 



dv[a(v) + ^(u)]Q' 1 (u - v) - 2Q 1 (u - t )a(t ) 

to 

r d 

[<j{u) + v(u)]Qi(0) + / dv — [a(v) + v(v)]Qi(v - u) 

Jto dv 



(42) 



since a(t ) = v{t$) = 1 . Now 



du[a(u) - v(u)}[a{u) + i/(u)]Qi(0) = (43) 



as cr 2 (u) — v 2 (u) = 0, for all u . Therefore the influence 
function reduces to 



F[a, u] = exp I du 

\J to Jto 



dvf(u,v] 



(44) 



with 

f(u,v) = -^[a^-viu^W^ + v^Qriv-u) 



4 
Ah 



a(u) - v(u)][a(v) - v{v)]Q 2 {u - u)(45) 



which may be compared to equation 7 of [T^l ] . 

We note (see also the appendix [Bjl that the second 
factor of (|45|) can be viewed as the effect of a classical 
fluctuating force. We will specialise later to the ohmic 
case when 



J{uj) « rjwe-"'" 



(46) 



for which Qi(u) will be well approximated by a delta 
function r]S(u). 

We can now follow Waxman ^| by differentiating the 
above propagator K to obtain a differential equation for 
the time evolution of K differing from his equation (15) 
only in the replacement of e by e(t) and by the replace- 
ment of his term in Q\ by one arising from the fixed spin 
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initial condition: 



iri—(aia 2 t\a 3 a i t ) 



where 



-^-[K(-aia 2 t\cr 3 a4i ) - K{a x - a 2 t\a 3 a4to)] 



he(t) 



((Ti - (T 2 )K{fj 1 a 2 t\a 3 cr i t ) 



(|) 2 (<ti - a 2 ) jf dwQxft; - t)^ < + > 
i (y) - cr 2 ) / dvQ 2 (v -t) < a(v) - > . 



(47) 



< tr(«) ± >= / d[<r] / dM\[(7]i;[i/]F[(T, v][<r(v) ± 

Jcr 3 ta Jcnto 



(48) 



r 



For the driven spin-boson system with factorizing initial 
density matrix and equilibrium fixed spin initial condi- 
tion the above integro-differential equation l|47l) is exact, 
and the first main result of this Paper. 



III. 



BLOCH-REDFIELD EQUATIONS FOR THE 
WEAK DAMPING CASE 



We now go on in this section ( IIII Afl to demonstrate 
that one may obtain a simpler set of equations-of Bloch- 
Rcdficld form-for the spin vector a by making a weak 
damping approximation, replacing terms under the re- 
tarded integrals of the propagator equation l|47|) by their 
equivalents for an uncoupled spin a^ ). We then ( IIII C|) 
identify the dissipation term x arising from the imagi- 
nary factor in the influence functional. In addition we 
identifvf llllDfl the terms due to the fluctuating force al- 
lowing us to complete the Bloch-Redficld equations, and 
compare them with those of Hartmann et al. 



< cr(v) ± v(v) >~< <y(v) ± v(v) >o This is in gen- 
eral a weak damping approximation. In the ohmic case 
it implies a <C 1 where a is the dimensionless fric- 
tion constant for the problem, a = riq^^irh and r\ 
is as defined in l|46|) . We note that linearity implies 

< a(v) ± v{v) >(o)=< °(v) >(o) ± < v{v) > (0) . 
Then < a(v) > =< (Ti|C/ (0) (t, v)a z U^{v, t )\a 3 >< 
<J 2 \U^ ' (t, tn)\<74 >* (see also equations (20) and (21) of 
|l5j\ Here U 1 -- ^ is the time evolution operator for a spin 
uncoupled to the bath: 



H s (t')dt' 



(49) 



A. Weak damping approximation 

We will now approximate the terms < cr(v) ± v{v) > 
in (|47|l by their value for an undamped system i.e. 



where T denotes the time-ordering operator. 



Using the definition of the inverse of we find 



< <t(v) >0=< 0"! 



\uW(t,v)a z uW-\t,v)UW(t,v)UW(v,to)\a 3 >< a 2 \U^{t,t )\a A >* (50) 
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but U^(t,v)U^(v,t ) = U(°\t,t ) so < we have: 
a(v) > =< a^^v^U^^vW^^to)^ >< 
(T2\\U( \t, to) |<X4 >* . Inserting a complete set of states 
of a z : 

1 = ^ |(J >< cr| = | + (Jl >< +£71 1 + | - £71 >< -£7l| 

(51) 

I 



< a{v) > = < o 1 \U^{t,v)a z U^- 1 {t,v)\a 1 > K <°) {a 1( r 2 t\a Z G A t Q ) 

+ < a 1 \U (0 \t,v)a z U^~ 1 (t,v)\-a 1 > {-a^t^a^) (52) 



where the superscript on K^°> indicates that it propa- 
gates the uncoupled spin density matrix. If 

tf(£7i£7 2 t|£7 3 £74to) = / d[a] d[v] A n {u] A* m [v]F[a, v\ 

(53) 

then 

if(°)(£7l£7 2 t|£73£74io) = / ' d[o] f ' d[v]A n [o]A* m [v] 

i.e. the case F[a, v\ = 1 (no environment). This is impor- 
tant because K propagates ^(o^a^to) to p{o\(J2t) (with 
t > t ) while will propagate the same p(<T3<74to) 

to a different final density matrix which we have called 
p(°\a\02t). Neglect of this point would lead to the ap- 
pearance of terms 0(a 2 ) in the final equation of motion. 
We note that the propagator J in equation (22) of 
should, in the notation of the present paper, have been 
written as J(°). 

We now want to evaluate the terms 

<a 1 \U ( ^{t,v)u z U^- 1 {t,v)\a 1 > . (55) 

To do this consider the density operator p^°'(t) which is 
the solution of 

M§-/ 0) (t) = lH S7 pW] (56) 
I 



< a(v) > = (fi(t,v) -if2(t,v)a 1 )K w {-a 1 <j 2 t\<T 3 a 4 t ) + £7i/ 3 (t, u^Vi^ks^to) (60) 

I 

This step is crucial because it allows the use of the spin into the equation for dtp{o\(J2t) will now give a soluble 
polarization vectors for an uncoupled spin. expression for the reduced density matrix /5(<7i a^t). Us- 

Substitution of our expression for dtK(tj\<j2k\ a 'i' J 'do) ing equations Q, i|60|) and 147( 1 we have: 



I 

Then pC°)(t) = U^(t,t )p {Q) {to)U^- 1 (t,t ) allowing us 
to form the identity 2p { -°\t) - 1 = U^> {t,t )[2pM (t ) - 
l]U^~ 1 (t,to). Defining a polarization vector a^°\t) by 
p (0) (t) = |(1 + a(°)(t) • a) we can write 2p - 1 as 

o<°> (t) • a = U<® (t, t ) [oM (t ) • ^C/W- 1 (t, t ). (57) 

If a(°^(to) ■ er = <j z which is the case we need in order to 
evaluate (|55|) . we can then (c.f. define an auxiliary 

polarization vector f(t,t ) by /(t, to) = a'°'(t). The 
dependence on to enters via 1(5 7(1, /(t, to) corresponds to 
the polarization vector for an isolated spin evolved from 
time to to time t subject to the initial condition that its 
polarization vector was (0, 0, 1) at time to- We will later 
use the equivalent notation w ^ (t, to) for this function, 
where the above initial condition is implied. 

Having established the meaning of the terms < 
<Ji\U(t 1 v)<T z U~ 1 {t 1 v)\(Ji > we are in a position to substi- 
tute the expressions for them into the equation of motion 
for K{a\a2t\<J-j,a 4 to). We first note that: 

<G 1 \U{t,v)a^\v)-aU- 1 (t,v)\a 1 >=£7 1 / 3 (t,«) (58) 

and similarly 

< £7i|C7(t, 1 ;)a( )(u)- CTt 7- 1 (t ) i;)|-£7i >= h{t,v)-ia 1 f 2 {t,v) 

(59) 

Hence equation (|52|l becomes: 
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iK—p(aia 2 t) 



= ih—^ 2J K{(J 1 G 2 t\a z a^)p{a z a4 a ) 



= ih 22 d t K(ai<7 2 t\(Jz(Jit ii )p{az<7ih) 



(?)< 



(Tl-(T2j 

x [ / dwQxCw - t)^-{/i(t,t;)p (0) (-Ti^2t) + fi(t, v)pW{a x - a 2 i)} 
J to dv 

dv Q 1 (v - t^i-ihit^icnp^i-awt) - a 2 p {0) (a 1 - a 2 t))} 

ft d 
+ (<?!+ ° 2 ) / dvQ^v-^—hit^Yp^Xa^t)] 
Jto dv 

-i [ dvQ 2 (t-v){fi(t,v)pW(-aia 2 t)-h(t,v)pW(a 1 -a2t)} 

Jto 

rt 

dvQ 2 {t - «){/ a (t, v)a 1 p^\-G 1 a 2 t) - f 2 (t,v)a 2 p^(a 1 - a 2 t)} 

'to 

'i{ai~<T 2 ) f dvQ 2 {t-v)f 3 (t,v)a 1 p ( - 0) {a 1 cr 2 t)}+O{a 2 ). 



(61) 



B. Identification of spin precession terms 



C. Identification of dissipation terms 



We note first that, after defining the reduced density 
operator p{t) and its corresponding polarization vector 
a(t) through 

p{cx 1 a 2 t) =< a x \p\a 2 >=< ai\hl+a(t)-a)\a2 > ■ (62) 

the terms in A and e(t) on the right of Q61JI will give the 
familiar spin precession term 



da 
~dt 



b(t) A a - 



We have also defined b(t) = (— A,0, e(i)) so that 



H S = -b(t) ■ a. 



(63) 



(64) 



We now consider the next three groups of terms, i.e. 
those in Q%, in equation (|61|l . These are due to dissi- 
pation and one may call them the systematic terms by 
analogy with the Fokker-Planck equation. The last is 
zero as (oi — cr 2 )(<Ti + 02) = o\ — o\ = for all a%, a 2 
(remember <7j = ±1 for i = 1, 2). The first can be found 
by forming commutators: 



(ai -a 2 )[p^{~a l a 2 t)+p^{a 1 ~a 2 t)] = < a x \[a z , a x p^(t)]\a 2 > 

+ < cri\[(J z ,p^(t)a x }\cr 2 > 



(65) 



r 



Now p^(t) — 5(1 + a^(t) ■ rr) so the above expres- term gives — — (J 2 )[<Jip(— a\G 2 t) — a 2 p(a\ — a 2 t)\ =< 
sion becomes < <j\ \2ia y \a 2 > . Similarly the second 
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(Ti| — 2ia x 1 02 >. Thus if 



^/i(M) = A(M) 



and similarly for and /g we have 



* 9 1 
— a ■ a 

dt2 



dv Qx(t - v)2i{f 1 (t, v)a y - f 2 (t, v)a x } 



to 



(66) 



where we have used the evenness of Q\{u) = Qi(—u). If 
e y is a unit vector in the y direction then <x • e y = a y etc., 

and defining a vector / (t, v) by 



/ (t,v) = (f 1 J 2 J 3 ) = —f(t,v) 



(67) 



D. Identification of fluctuating force terms 



We finally have the group of terms in Q 2 from equation 
(ED- First 



we then have 



da 
~dt 



b(t) Aa + X + 



with 
X 



(|) / dvQ 1 (t-v)e z Af'(t,v). 

J to 



(68) 



(69) 



- < (Ji\a y 0) (t)[a z ,[a x , a y ]}\a 2 > 

+- <<Ti\a®\t)[cr z ,[<T x ,cr z ]]\<T2 > 
0+ < o^^ifa^ > 



(70) 



This is completes the derivation of the dissipative part 
of the Bloch-Redfield equations for this initial condition. and similarly 



J 



i(ai - a 2 )[<Jip {0) (~<7 1 (T2t) + (T2p {0, {ai - a 2 t)] =< (Ti| - 2a y a ( " ) (t)\a 2 > 



(0), 



(71) 



while 

{(Ji-a 2 ) 2 p{(Ji<J 2 t) =< a 1 \{2a x a^\t)+2a y a y 0) (t)\a 2 > . 

(72) 

If we substitute these values into equation iJSTJl and use 
the auxiliary vector / we find that the equation of motion 
for a(t) is: 



+ J(f ) 2 dvQ.it -v)e z A f'(t,v) 

+Gz A \ (f ) 2 / dV Qa {V ~ <)/(< ' V) A a(0) {t) 
= bAa + x-e z A{tl>{t)Aa {a) ) + 0{a 2 ) (73) 



where we have defined the retarded integral 

*j = -\ (f Y j[ dv Q 2 {t - v)Mt, v) (74) 

Equation (|73|) has been derived by omitting terms of 
0(a 2 ) and higher. While this procedure seems reason- 
able it leads, at long times, to unreasonable results in, 
for example, problems where b is independent of time. 
In these cases the term continues to oscillate in- 
definitely, and thus the system will not tend to a time 
independent state of equilibrium, as we believe it should. 
The simple replacement of a' ' by a appears to cure this 
problem; for a static b the system does equilibrate [l5| . 
We view this as a selective resummation (to infinite or- 
der in a) that is sufficient to guarantee more correct long 
time behaviour, which we make due to its physical rea- 
sonableness. 
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We thus replace a°(t) in equation l|73[) by a(t) to give: 
3cl 

— = b(t)Aa + X -e z A (-0(0 A a) (75) 
at 

Equation l|75|) and its derivation by path integral, with 
the auxiliary definitions of x ( I69|) and ip ( I74|l forms the 
second main result of our work. 



Comparison with Bloch-Redfleld equations of 
Hartmann et al. 



Hartmann et aL[lJ] quote a set of Bloch-Redfield equa- 
tions obtained from projection operator methods as their 
equation (4): 



O ti 



+ 



e 




\ 




-e 




A 


1(5. 


- 


A 


) 




1 XX 







-T xz \ 







r 












o / 




However the difference would disappear if we were to 
revert to a factorizing initial condition without fixed spin. 
The other point to stress is that both our derivation and 
theirs arrive at a fluctuating force term with a , the sub- 
stitution of a at that point has to be an additional as- 
sumption. 



IV. BLOCH-REDFIELD EQUATIONS IN THE 
OHMIC AND ZERO TEMPERATURE CASES 

We now use the result of Waxman and Zhang [l^ 
that in the ohmic case the full expression for the time 
dependent x vector ( EHl greatly simplifies to the time 
independent x^x- In the critical damping case studied by 
Shytov y| where J(u>) ~ w _1 both x and y terms would 
be present in X- 



A. Simplification of the dissipation term for ohmic 
spectral function 

In the ohmic case, equation (|35[l becomes Qi(u) w 
rj5{u) so 



Ay 





(76) 



x = fkA/'(M) 



(80) 



where <J x .y,z are spin expectation values i.e. a x , y . z in 
our notation. To establish the correspondence with our 
work note first that they have a complex correlation func- 
tion: 



where the extra factor of 1/2 comes from the use of the 
delta function at an endpoint of the integral. Now 

f'(t,t)=]im^-f(t',t)=bAf(t,t) (81) 



M = - 

7T 



where the last step follows from the fact that 



dwJ(w)[cosh(?iw/2fc B T 



= M + iM 



iut)/ smh(huj/2k B T)} 
(77Xs 



which is equivalent to Qi + iQ[ in our notation. They 
also define transition amplitudes for an uncoupled spin 
between right and right, and left and right, pairs 
of states as U RR =< R\U(t,t')\R > and U RL =< 
R\U(t,t \L > respectively, and an auxiliary function 
F(t) = 2 J" * dt ' M" (t - t')U RR U RL so that A x = ImF(t) 
and A y — ReF(t). Finally they have the damping tensors 



— f(t',t)=bAf(t',t). 



lim /(*',*) = /(*,*) = (0,0,1) 



we have 

b A /(*, t) = (-A, 0, e(t)) A (0, 0, 1) - Ae y 



(82) 



(83) 



(84) 



so 



e z A (6A/(t,t)) = Ae z Ae % 



Ae, 



(85) 



Fij defined by I\y = f Q dt' M' (t-t')bij(t, t). If we follow We thus find that the dissipation term ||HU|) becomes: 



Hartmann et al. by defining b yz = —2Im(U RR U R l) and 
b xz = —2Re(U RR U R l) , it becomes clear that their b yz 
corresponds to our a° and b xz to our a° x . 

We thus find that their T xx — T yy = —ip z , and also 
T xz = ip x , T yz = ip y . This gives an identical structure to 
ours for the fluctuating force term, while their dissipation 
term A — A x e x + A y e y (= x) is different, being 



ml 

2% 
—naAe. 



("A)e a 



A x = 



A, 



dt Q[(t ~ t )a y 



la 



dtQ l (t~t)al 



(78) 



(79) 



= Xe x (86) 

We note that this term acts in the negative e x di- 
rection. By taking the dot product of a with (|75Jl to 
find da 2 /dt we see that the dissipation term thus tends 
to reduce the length of the vector a, in contrast to the 
magnitude-conserving spin precession term. This is di- 
rectly analogous to the way that an electric field can 
change a speed while a magnetic field cannot, in the 
Lorenz force law. 
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B. Additional simplification of fluctuation term in 
zero temperature limit: 



A. Illustrative numerical solutions of ohmic 
Bloch-Redfield equation at zero temperature 



We have, provided that ui c t 3> 1, 



Q2{V) = -\-J-Zl 



71" (W~c 



2\2 



i._(l sinh ^)- 2] (87) 



This has two components corresponding to zero and 
finite temperatures respectively. We will here specialize 
to the zero temperature limit. In this case ^3 becomes 
— 2ira while ip x and tpy tend to zero. 



V. APPLICATION TO 
LANDAU-ZENER-STUCKELBERG PROBLEM: 



The dissipative LZS problem as usually considered 
(e.g. 01) corresponds to the driven spin-boson Hamil- 
tonian with the specific choice e(t) = et, although other 
choices may be of interest (e.g. [23). We emphasise that 
in this paper we have kept general time dependence e(t) 
of Q up to now. The standard LZS choice corresponds 
to a spin modelling a particle in a biased double poten- 
tial well for which an external driving force acts to reverse 
the sign of the offset between the wells at t = 0; t hav- 
ing increased from a large negative starting value t = t$. 
The quantitities of interest are usually the long time, e.g. 
t — > 00, values of the probability P^z that the spin is in 
the "down" (—1) orientation if started "up" (the conven- 
tional "tunnelling probability") or, equivalently, the long 
positive time expectation values a x ,a y , a z of the spin op- 
erator a. 

In absence of coupling to environment the LZS tun- 
nelling probability can be obtained exactly as (e.g. 



Pi 



LZS 



1 



-7rA 2 /2e _ y 



where t z — A/e is the Zener (characteristic tunnelling) 
time. 

Plzs is ~ 1— e~ s / h where the action S is the energy for 
one hop times the Zener time, which, to within a numer- 
ical factor, can be derived from an instanton calculation 
0. The inverse of the tunnelling angular frequency ex- 
pressed in units of r z is 2/(Ar z ) = A, a nonadiabaticity 
parameter [2^, which increases as the speed of change 
of the energy levels increases i.e. as the system becomes 
more nonadiabatic. 

Figures 1 and 2 illustrate the Landau-Zener behaviour 
(see also [13 )• For large A the value of S tends to zero, 
as does Plzs and so the probability of not tunnelling 
tends to 1 i.e. in the fast passage limit the system does 
not tunnel. In the opposite limit of adiabatically slow 
change, the particle clings to the instantaneous eigenstate 
and must tunnel, so Plzs tends to 1, i.e. the probability 
of not tunnelling tends to 0. 



To solve our Bloch-Redfield equations we need to 
rescale time to the Zener time t z . After doing this the 
ohmic version becomes: 

^ = BAo(T) + y ei -e 2 AWT)Aa(T)) (89) 

We first note that this form shows why the effect of the 
dissipation term is negligible for the (fast passage) non- 
adiabatic limit of large A. It is simply because the dis- 
sipation term is inversely proportional to A. This means 
that (within the limits of our approximation) the faster 
a damped qubit is switched, the less it will be affected. 

We show representative solutions of (|89l) in the ab- 
sence of fluctuations as an illustration of the size of the 
perturbation of the probability caused by the spin-bath 
coupling. Figure 3 shows how (1/2)(1 + a z ) changes in 
the most adiabatic case studied, A = 1, comparing the 
uncoupled evolution (blue curve) with a = 0.001 Figure 
4 is for the same case as Figure 3 but shows the 3 com- 
ponents of spin, only for the damped case. Figures 5 and 
6 are analogous but show the fast passage A = 15 case. 
Figures 7 and 8 are also for the fast passage case A = 15 
but are for stronger coupling to the bath a = 0.01. At 
this level of coupling the damping of quantum oscilla- 
tions is very noticeable, although it does not alter the 
final tunnelling probability. 



VI. CONCLUSIONS 



This paper has extended the static spin boson study of 
Waxman [l5| to the driven spin boson case, obtaining an 
exact integro-differential equation for the time evolution 
of the propagator of the reduced spin density matrix, the 
first main result. By specializing to weak damping we 
then obtained the next result, a set of Bloch-Redfield 
equations for the equilibrium fixed spin initial condition. 
Finally we showed that these equations can be used to 
solve the classic dissipative Landau-Zener problem and 
illustrated these solutions for the weak damping case. 
The effect of dissipation was seen to be minimised as 
passage speed increased, implying that qubits need to be 
switched as fast as possible. 
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{x a } (or x) set of harmonic oscillator positions x a 
Zq partition function 

z\ , z 2 temporary sum and difference co-ordinates 



APPENDIX A: NOTATION 

a. dimensionless coupling constant 
a(t,to) spin polarisation vector 

a,(°\t, to) (= f(t, fo))spin polarisation vector in uncou- 
pled case 

A[a] Spin amplitude functional for path containing n 
flips 

A[(i,f] Spin amplitude in stochastic case 

Inverse temperature 

b vector derived from 2-level Hamiltonian 

{c a } Set of spin-bath coupling constants c a 

d[x] measure for path integral 

d[a] measure for spin path integral 

A tunnelling matrix element 

D Correlation function in stochastic case 

e(t) bias 

?i ohmic coupling strength 

/ Fluctuating force 

F[a, v\ Influence functional 

Fa=o Influence functional in uncoupled case 

e x ,e v ,e z unit vectors 

h Planck constant 

H(t) Hamiltonian with corresponding action S[a, x] 
H a (= H(t = to)) 

HB{t) Bath Hamiltonian with corresponding action 

S B [x] 

Hj(t) Interaction Hamiltonian with corresponding 
action Si [a, x] 

Hs(t) Spin Hamiltonian with corresponding action Ss[<j] 

J(u>) bath spectral function 

K(°i (t) spin propagator for 

K(t) spin propagator for p 

X dissipation term 

ks Boltzmann constant 

A adiabaticity parameter 

m a oscillator mass 

u) a oscillator frequency 

lo c bath cutoff 

ipi i = x,y, z 

Plzs Landau-Zener tunnelling probability 

P[f] stochastic force probability functional 

qo inter-well spacing 

Qi retardation function 

Qi fluctuating force correlation function 

p osc Single oscillator density matrix 

p density matrix 

p spin density matrix 

u x,y,z Pauli spin matrices 

t z Zener time 

t time 

to fiducial time 
T temperature 

U unitary time evolution operator 



APPENDIX B: IDENTIFICATION OF 
FLUCTUATING FORCE AND DERIVATION OF 
STOCHASTIC EQUATION OF MOTION FOR 
PROPAGATOR OF DRIVEN SPIN-BOSON 
SYSTEM 

We now show that the second term in the influence 
functional is equivalent to a random fluctuating force 
leading to an alternative equation of motion for the prop- 
agator of stochastic (quantum Fokker-Planck) type. This 
result is still quite general, but we confirm that the known 
delta-correlated force would be recovered in the high tem- 
perature limit for an ohmic choice of spectral density. 

We first define a probability functional P[f] which 
gives the probability that a member of a statistical en- 
semble experiences a force / at a given instant: 

P[f] = cxp-1 J ~f{u)D-\u,v)~f{v) (Bl) 

normalised such that 

d[f]P[f] =< 1 >= 1. (B2) 

Then we define an average over the ensemble of forces 

< e~ f f( u ^ h( - u } >J= I h(u)D(u,v)h(v) j-gg-j 

which implies that the autocorrelation function is given 

by 

</>)/>)>= D(u,v). (B4) 
If we expand D(u, v) in terms of its eigenvalues {A„} 

D ( u , v ) = ^n(j)n(u)<j) n (v) (B5) 

71 

we can define its inverse D^^.v) by 

D-\u,v) = Y J K 1 4>n{u)4> n {v) (B6) 

n 

where we require A„ > for all n. We define the corre- 
lation function Q2{u,v) through 



D(u,v) = 2hQ 2 (u - v). 
Hence if we take the average 



(B7) 



J* du[a(u)-v(u)]f{n) 

< e Jt o > f — 



-K^f J* du J' dv[a(u)-a(u)][a{v)-a(v)}Q 2 {u-v) 
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we obtain an object which has the the same form as the 
second, real factor in the influence functional (|45l) . The 
identification of the correlation function Q2 (it — v) can be 
made with our previously defined Q2{u — v) (|38|l provided 
that the Fourier transform is positive. We recall that 
Q2 is of the form of a cosine transform so we can see 
positivity is satisfied. 

We can thus define a new amplitude A[a, f] and an 
influence functional F\a, v\ via 



A[a,/] = (^rexp-i 



'0 



M«) (B9) 



and 



F[a,v] = ex P-r(f) 2 I du I ^[<r(u)-^(w)][a(u)+^(u)]Q / 1 (w-«) 

•J to J to 

(BIO) 

The stochastic analogue of (|47|l is, before averaging with 
respect to /: 



ihd t K(ai<72t\v3Vi) = 



'-^■[^(^ Cr l Cr 2i|c r 3CT40) - A'((Ti - cr2i|cr 3 <7 4 0)] 

+ + § /(*)) (ff X - a 2 )K(a 1 a 2 t\a 3 a i 0) 

+ (y ) 2 (<ri - a 2 ) J dvQ[(t -v)< a(v) + v{v) >, 



(Bll) 



where we denoted the stochastic propagator by K. This 
propagator equation could now be used with Q to obtain 
a quantum Fokker-Planck equation for the driven spin- 
boson model, analogously with 0, [2(J . 

The replacement of the Q2 term has thus contributed 
a fluctuating bias to the existing time dependent bias. 
In consequence the stochastic analogue of the Bloch- 
Redfield equation J7SJ), before averaging with respect to 
/, becomes: 



da ... 
w =b(t)A«- 



f f(t)e z A a 



(B12) 



This is a weak damping quantum Langevin equation. 
Here < f{u)f(v) > f— 2hQ2(u — v) and < a >j is the 

physical quantity. f(t) has the physically appealing in- 
terpretation of being the fluctuating force on the spin as- 
sociated with a given member of the ensemble of systems 
considered. We believe that this may be more tractable 
than the Bloch-Redfield form in cases when the retarded 
integrals would otherwise need to be evaluated. 

We note that we can recover the known delta- 
correlated behaviour of Q2 for the high temperature 
(t — > 0), ohmic spectral density (Eq. I46fl limit. Using 
(021 we nrs t substitute for J(w) in C I55) l. and then, pro- 
vided the bath cutoff frequency w c is sufficiently large 
but lot = TiLa/hsT still sufficiently small, we can approx- 
imate cothx by l/x. This gives us Q2 — [ 2t ^ t ](5(u) and 
thus 



< /»/>') >= 



A V k B T 



8{t - 1') 



(B13) 



APPENDIX C: GREEN'S FUNCTION SOLUTION 
OF OHMIC QUANTUM LANGEVIN EQUATION 

We here solve the quantum Langevin equation pertur- 
batively using the Green function method of Zhang |17| . 
We find the intriguing result in the ohmic weak damping 
case that the fluctuations enter only at second order in a, 
and so give an expression for the tunnelling probability 
to first order in a. 

We use the methods of equations (46-52) of ^3] 7 but 
unlike Zhang we do not differentiate the Green's function 
solution. 

We have: 



^=b(t)Aa + x+f/(*)e,Aa (CI) 

where /(f) is the fluctuating force. This allows us to 
write a solution of the stochastic equation as 



a(i,/) = M(i,t )a(i )+ / M(t,t') X (t')dt' (C2) 

(c.f. equation (46) of 01) w ^fi a Green function for 
the homogeneous equation 



^-bAa-|/(f)e 2 Aa = (C3) 

We then take out the evolution operator u(t, t') for 
the unperturbed spin, and then (equation (49) of [l7| 
rearranged) note that 



u(t,t')Af{t,t') = M(t,t') 



(C4) 
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giving 

a(i, /) = u(t, t )M{t, t )a(t ) + / u(t, t')M(t, t') x (t')dt' 

(C5) 

Now, following equation (52) of the formal solu- 
tion for Af is 

JV(t, t') = 1 + J / dt n F{t u t')F{t n ,t') 

(C6) 

I 



Now < Af(t,to) > is 

I 



to second order in J 7 (which we know will be first order 
in a). We may take < !F(t\,t') > to be zero for the 

I 



because [e z A (e z A]a(i ) = for the case we chose, where 
a(to) was (0,0,1). The terms under the Green's func- 
tion integral (i.e. eauationf IC8|) for N(t,t')x{t') are all 
second order in a except the first one which from above 
expansion can be seen to be just xC^Oi so we are with 
a very simple solution 

a(t) = u(t,t )a(t ) + [ u(t,t')x(t')dt' (C12) 
= a°(t,t )+ x / a 1 (°)(t,0* / (C13) 

J to 

where a 1 ' ' (i, t') is just the undamped solution a at time 



with 



T = ff(t)e z A . (C7) 



If we now average over the stochastic force, we will have 



(C8) 



(C9) 

I 

fluctuating force without serious loss of generality. 
We note that 



(CIO) 



(Cll) 



t, but started at (1,0,0) at t'. In the ohmic case after 
rescaling times by t z we have x = 27ra/A (c.f. 1551 1 so we 
have a linear dependence on a for a given A. Checking 
this prediction with numerical solutions for a = 0.005 
agrees very well with the above result. We note this is 
a consequence of the spin-up initial condition, as others 
might give a first order contribution from the fluctuating 
force. It is however independent of the spectral form 
of the environment once that initial condition has been 
specified. 



i(t) =< a(t, /) > f= u(t, t ) < Af(t, t ) > a(t ) + / u(t, t') < Af(t, t') > X (t')dt' 



to 



l+\ dh< F{ti,t') > 



dhdt 2 < T(t u t')T{t 2 ,t') > 



tQ JtQ 



<Af(t,to) > a(i ) 



[1+ / / dtxdt 2 <^(ii,i')^'(*2,t / ) >]a(*o) 

J t J t 

[ 1 + tI / / *i*a < /(*0/(*a) > e, A (g,A]a(to) 
n J to Jto 

a(t ) 
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FIG. 1: Dependence of the Landau-Zener-Stuckelberg tunnelling probability Plzs and the probability of not tunnelling 1 — Plzs 
on the adiabaticity parameter A in the absence of spin-bath coupling. 




FIG. 2: Evolution from t = — 10r z of 1 — Plzs to t = 10r 2 for spins started with a = (0,0, 1) in the absence of spin-bath 
coupling (after Berry |25||1. A increases as the curves go up the plot, taking values of 1, 2.5, 5, 15 
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FIG. 3: The evolution of (1/2)(1 + a z ) with time in medium passage case, A = 1, comparing the uncoupled evolution (blue 
curve) with a = 0.001 (red curve). 




FIG. 4: Three components of spin for the same case as Figure 3 showing only the damped cases. 




FIG. 5: As Figure 3 but for the fast passage case A = 15. 




FIG. 6: As Figure 4 but for the fast passage case A = 15. 
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FIG. 7: As Figure 5 but for stronger coupling to the bath a = 0.01. At this level of coupling the damping of quantum oscillations 
is very noticeable, although it does not alter the final tunnelling probability. 




FIG. 8: As Figure 6 but for stronger coupling to the bath a — 0.01. 



